CHERENKOV GLUE IN OPAQUE NUCLEAR MEDIUM 
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. Abstract 

(— I I The spectrum of Cherenkov gluons is calculated within the frame- 

Qh' work of in-medium QCD. It is compared with experimental data on 

p • the double- humped structure around the away-side jet obtained in the 

0) ; experiments STAR and PHENIX at RHIC. The values of the real and 

imaginary parts of the permittivity of the medium created in ultra- 
relativistic nuclear collisions are obtained from these fits. It is shown 
that accounting for an additional smearing of transverse momenta of 



> 



I Cherenkov gluons in the medium allows to provide a very good de- 

' scription of the experimental data. 
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Remarkable experimental results revealing an existence of double-humped 
structure of away-side azimuthal correlations in central high energy nuclear 
collisions obtained in the experiments STAR [H [2] and PHENIX P H |5] at 
RHIC have triggered a new wave of interest to a phenomenon of Cherenkov 
radiation of gluons first considered in [6l [3 E]. A simple model of scalar 
excitations in strongly interacting quark-gluon plasma leading to Cherenkov- 
like pattern was considered in [H] . An alternative explanation of the observed 
double - humped structure in terms of Mach cones was discussed in [TOl 
[TT] . For a recent comparison of these two mechanisms see [12]. Let us 
note however that a recent detailed study [13] presents strong arguments 
against the Mach-cone interpretation by proving that a dominating feature 
of hydrodynamical picture of energy loss is a strong diffusion wake leading 
to a pronounced maximum in the forward direction. 

In an ideal case of transparent medium Cherenkov gluons are radiated at 
fixed angle with respect to the direction of propagation of the color current. 
Their spectrum is described by the well-known Tamm- Frank formula [131, [15] . 
The experimentally observed double-humped angular distribution shows the 
structure that is noticeably spread around the angles corresponding to the 
maxima. Among the physical reasons causing such spreading are, in particu- 
lar, the opacity of the color medium and the finite propagation length of the 
color current. Some effects were treated in [T6j. In the present paper we con- 
sider a model combining an account for the opacity of the color medium and 
an additional smearing of "pure" Cherenkov pattern in the plane orthogonal 
to a propagation direction of the radiating current that allows to achieve a 
good description of the available experimental data. 

The equations of in-medium QCD were proposed in [17]. At tree level and 
in quasi-abelian approximation they coincide (up to color factors) with the 
corresponding equations in QED. Therefore, in this approximation Cherenkov 
gluons are described in complete analogy with Cherenkov photons and cor- 
respond to the specific singularity of an abelian Lienard-Wiechert potential 
having its origin in the nontrivial (chromo) dielectric permittivity e. In [17] 
the simplest case of color singlet constant real permittivity was considered. 
In the present paper, still working in the quasi-abelian approximation and 
considering the color singlet permittivity, we shall consider a more general 
case of opaque medium characterized by complex permittivity 



where ei{uj) = ei(— cj), €2(0;) = — e2(— u;), uj is the gluon energy. For simplic- 
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ity we neglect the spatial dispersion. 

Truly non-abelian effects can be considered by taking into account the 
nontrivial tensor chromoelectric permittivity in the color space and the non- 
linear terms in the Yang-Mills equations in the medium [T7j. These effects 
will be considered in the future publication [18]. 

To develop a description of the angular pattern observed in the flow of 
final hadrons in high energy nuclear collisions [H El El IH IS] in terms of the 
Cherenkov gluon radiation one has to 

(a) Describe a kinematical pattern characterizing the initially produced 
hard partons serving as colored sources of gluon Cherenkov radiation taking 
into account the experimental cuts on pseudorapidity and transverse momen- 
tum PE]. 

(b) Write down a spectrum of Cherenkov gluon radiation in the opaque 
medium for the above-described color sources. 

(c) Describe a conversion of Cherenkov gluons into observed hadrons 
taking into account the experimental cuts on transverse momenta of final 
hadrons [21 E] • 

These three steps are realized in a Monte-Carlo procedure. In total the 
model includes three parameters, described in more details below - two re- 
lated to the Cherenkov gluon radiation, the magnitudes of the real ei and 
imaginary t2 parts of the permittivity of the medium, and a scale Aj^ con- 
trolling the transverse smearing of Cherenkov gluons . The values of these 
parameters are fitted to achieve the best possible agreement with the exper- 
imental spectra [21 E]- 

The description of the initial hard color sources was performed through a 
Monte Carlo simulation of hard parton-parton scattering at RHIC energies 
with PYTHIA [TH]. An initial pool of two-parton configurations resulting 
from these scatterings consisted of those in which the trigger jet gave rise to a 
pion hitting the trigger intervals in rapidity, | rj\ gjig^j, < 0.7 and | rj\ phenix < 0.35, 
and transverse momentum, where we have chosen the intervals 3 GeV < 
I pi^^l < 4 GeV for STAR [2] and 2 GeV < | p|[^| < 3 GeV for PHENIX [5]. 
The fragmentation of gluon (quark) generating the trigger near-side jet into 
pions was described by standard fragmentation functions measured at LEP 
[20j . For simplicity we have considered only the dominating subset of events 
with gluonic away-side jet. 

Let us now turn to the differential Cherenkov energy loss spectrum in 
the opaque medium per length dl along the direction of the radiating gluons 
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(quarks) [211 [22] 



dW _4asCv(F) cose{l-cos^e)T{uj) 

(2) 



Lodldudcpd cos 9 vr (cos^ ^ - CM)^ + ^(a;) ' 

where Cv{f) are the Casimir operators for the adjoint and fundamental rep- 
resentations of the color group correspondingly, 6, are polar and azimuthal 
angles with respect to the emitter propagation direction and 



(3) 



so that C{^) controls the location of the maximum and T{uj) controls the 
spreading around it. This a'la Breit-Wigner shape replaces the 5-functional 
angular dependence characteristic for real e. A derivation of Eq. ([2]) in terms 
of the energy flow of radiated gluons is given in Appendix 1. For (£2/61)^ <^ 1, 
the energy-dependent maximum of the differential spectrum ([2]) is at 



The second relation is just the Tamm- Frank expression p3l [15] for ultrarel- 
ativistic particles. Let us stress that in Eq. ([2]) the angle 6 of emission of 
Cherenkov gluons is measured with respect to the direction of propagation of 
the radiating color current which in the considered case is that of an initial 
gluon of the away-side jet. In the ensemble of initial jets the near-side and 
away-side jets are in general not angularly balanced in the laboratory system. 
However, it was argued in [16] that an angularly balanced configuration is 
in fact dominant because a generic shape of structure functions leads to a 
suppression of energy-mismatched configurations of the initial partons and 
consequently of the angular mismatched configurations of the final ones. In 
our Monte-Carlo program this mismatch is, nevertheless, taken into account 
by using PYTHIA. 

At the second step of our Monte-Carlo procedure we generate Cherenkov 
gluon radiation by the away-side gluon jet. To illustrate the correspond- 
ing kinematics let us consider a simple example of the Cherenkov spectrum 
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generated by the source having a direction of propagation orthogonal to the 
colhsion axis z that allows an analytical treatment. In order to compare the 
spectrum ([2]) with the experimental one in this kinematics we should rewrite 
it in terms of experimental laboratory polar and azimuthal angles Qj, and 
It is easy to see that for the above geometry 

cos6' = I sin6'2, cos0l|, 

for (^i^ counted from the away-side jet. 

The number dN of emitted gluons per length dl with energy and angles 
within intervals dio dO^ d(j)L is 



dN 



dl duj d(f)L d cosOl 



AasCv I sinOi cos0l|(1 — sin^ 9l cos^ (f>L)T{u!) 
(sin^^iCosVi -CM)Vr2(cu) 



(5) 



The experimentally observed hadron spectrum results from hadronization of 
the gluon distribution dN/dcpL obtained from Eq. ([5]) by integrating over 6l: 



dN 



du dl ( 
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A{uJ,(f)L) 

B(u;,(Pl) 



T(uj) cos^ < 



(7) 



The gluon spectrum (jSj) reveals the double- humped structure of interest. 

Let us stress once again that in our Monte-Carlo procedure we consider a 
more general situation where the direction of the gluon generating the away- 
side jet is fixed within each configuration satisfying the above-described trig- 
ger conditions imposed on the properties of the near-side jet in experiment. 

To realize a Monte-Carlo procedure of generating the Cherenkov spectrum 
we have to specify the functions ei{uj) and 62 (c^;). Taking into account that a 
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range of (laboratory frame) transverse momenta of the away-side hadrons in 
which the double-humped structure is observed is quite narrow ( ~ 2 GeV) 
and the binning of experimental data we aim to reproduce was performed at 
a scale of 1 GeV we can with a good accuracy neglect the effects of dispersion 
and consider energy-independent ei^2 = const. The values of these constants 
are determined through fitting the experimental data. Within this assump- 
tion the spectrum of produced Cherenkov glue is simply energy-independent, 
dN/dujdl = const. Realistically the possibility of Cherenkov emission is re- 
stricted to some finite interval of energies u < cjmax so that 

dN ^, , , , 

OC (cUmax " t^j , (8) 



duj dl 



where oApax is the highest characteristic resonance excitation energy of the 
mediuno. In our simulations we have chosen cUmax = 3.5 GeV. We have 
verified that our results are in fact weakly sensitive to the exact value of 

The thus generated flow of Cherenkov glue should, of course, be trans- 
formed into that of final hadrons. There exist several phenomenological 
schemes describing this conversion. In our case it is convenient to use a lan- 
guage of fragmentation functions Dg{x, p± \ Q"^) which generically describe a 
probability for a gluon with energy E and virtuality to produce a hadron 
with the energy xE and transverse momentum measured with respect to 
the direction of propagation of the initial gluon. 

One is first tempted, relying on the soft blanching hypothesis, to ascribe 
the same shape to the spectrum of hadrons resulting from fragmentation of 
the Cherenkov hadrons which amounts to assuming Dg{x,p± \ Q"^) oc (5(1 — 
x) as supported by the experimental evidence obtained in the physics of 
multiparticle production from e~^e~ — > jets. However, this does not lead to 
the fully satisfactory description of experimental data because of predicting, 
in contradiction with the experimental data, that the probability of radiating 
Cherenkov glue at angles | 0l| > n/2 with respect to the away-side jet is 
strictly zero. These "dead zones" appear because the Cherenkov gluon "halo" 
can not spread to angles larger than 7r/2 to the direction of the emitter. It 
is clear that such "dead zones" will in fact be present for any fragmentation 
mechanism that does not generate transverse momentum with respect to 



^For a discussion of interrelation between Cherenkov gluons and hadronic resonances 
in the medium see 1^. 
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the direction of propagation of the initial Cherenkov gluon which in our 
case corresponds to considering a p^-independent fragmentation function 
Dg{x, I Q^). We confine our consideration to fragmentation to hght hadrons 
because experimentally the share of other bosonic resonances is quite small 
(in particular, p : oj : 0=10:1:2 according to [2l])- In our analysis we used a 
simple parametrization of the fragmentation function (x, px| Q^) including 
the effects of transverse smearing: 

/^J(a;,pJQ^)oc(l-xf-^=i=exp|-^} , (9) 

where the x-dependence is that of a gluon fragmentation function at the 
reference scale Qq = 2 GeV [20j . 

Thus, we have a set of three parameters: ei, t2 and the transverse smearing 
scale A_L. 

The values of the parameters that provide the best description of exper- 
imental data of STAR and PHENIX collaborations for the ranges of trans- 
verse momenta in the associated jet 1 GeV < | p^^*^! < 2 GeV for STAR and 
2 GeV < I pI^'^I < 3 GeV for PHENIX are shown in Table 1. 



Table 1 



Experiment 


^max 


ei 


62 


Ax 


STAR 


1.04 rad 


5.4 


0.7 


0.7 GeV 


PHENIX 


1.27 rad 


9.0 


2.0 


1.1 GeV 



Let us stress that a difference in the fitted values of ei for STAR and 
PHENIX originates from different positions of angular maxima 6'max in the 
corresponding experimental data. (May be, it could indicate on necessity to 
take the dispersion into account.) However these values are quite stable in 
both above approaches because they are defined by maxima positions but 
not by the humps widths. 

The value of €2, on the contrary, is influenced by the widths. Nevertheless, 
the ratio (£2/61)^ < 0.05 stays quite small in all cases. 

As to the transverse smearing parametrized by A^, one can get a feeling 
of the magnitude of this effect by considering a typical distance covered by 
the away-side parton L ^ 5 fm [25j . This gives a smearing rate per fermi of 
Ai//m ~ 0.2 [GeVVfm] . 
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A^=0,7 GeV/c 



Figure 1: Away-side azimuthal correlations for STAR, circles - experimental 
data, triangles - simulation. 



The resulting angular spectra for STAR and PHENIX are shown in Figs. 
1 and 2 correspondingly. 

We see that the positions of the maxima (and therefore the values of 
ei) are quite stable to accounting for additional smearing on top of that in 
Eq. ([6]). However, it is indeed important in achieving a good description of 
the widths of humps in experimental data as seen in Figs. 1 and 2. The 
shape of humps in the former "dead zone" determines the parameter Aj^, 
which, in its turn, influences 62- 

To conclude, the values of the real part of the nuclear permittivity ei found 
from the fit to experimental data of RHIC (albeit somewhat different in the 
two experimental sets) are determined with good precision. Their common 
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A 
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Figure 2: Away-side azimuthal correlations for PHENIX, circles - experimen- 
tal data, triangles - simulation. 

feature is that they are noticeably larger than 1. This shows that the density 
of scattering centers is quite large and the nuclear medium reminds a liquid 
rather than a gas (for more details see [IS])- The accuracy of the estimate 
of the values of €2 is much less but more important is the fact that they are 
rather small compared with ei. 

In this paper we have described a qualitative model that allows, in a 
framework of a reasonable phenomenological picture, to reproduce the exper- 
imentally observed double-humped structure of the away-side jets in heavy 
ion collisions at RHIC [H [21 [3l IH [5] . To develop a quantitative description 
one needs to consider a more specific microscopic model for chromodielec- 
tric permittivity and study the effects of hadronization in more detail using 
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Monte-Carlo simulations. Such studies aimed at working out a quantitative 
explanation of the dependence of the effect on the range of triggered trans- 
verse momenta, in particular a disappearance of the effect at increasing j9_l, 
and a detailed description of the central domain related to fragmentation of 
the radiating particle are now in progress. Let us stress that in working out 
such quantitative description it is crucially important to have experimental 
data in a broader range of (laboratory) transverse momenta, in particular on 
configurations with a pronounced mismatch between the momenta of trigger 
and associated jets. 
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Appendix 1 

The energy losses due to emission of transverse gluons read 
dW iCvg^ [ , f „ uj - (kv)2 



dz 2n^ 



where v is a unit vector determining the direction of propagation of the radi- 
ating color current. In cylindrical coordinates {k, q, 0) with v = {v^, Vy, v^) = 
(0, 0, 1) we have 

dW iCvg^ /" J J 3 P J, 1 1 

(11) 

An angle at which gluons are emitted can be found from the ratio of the 
corresponding components of the Poynting vector 



SI El{u,\^) k^- 



^e.^.^.^. (12) 
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A short calculation leads to the following expression for the differential spec- 
trum of emitted gluons 

where 

p[tan^|CM,rH] = ^ T-^TT^^T^T-^iTT (14) 

and 

C(a;) = ei(cj) - 1, 

fH = esH- (15) 

Finally, a transformation to the variable cos^ 6 leads to the spectrum of 
Eq. © [an [22] 

1 dV^ _ 2Cy^^ cosg(l -cos^g)rH 
udzdudcose~ vr (cos2 ^ - C(cu))^ + r2(cu) ' 
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